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ABSTRACT 

Knots are commonly found in nearby planetary nebulae (PNe) and star forming 
regions. Within PNe, knots are often found to be associated with the brightest parts of 
the nebulae and understanding the physics involved in knots may reveal the processes 
dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an 
ideal target to study such small-scale (^300 AU) structures. We have obtained infrared 
integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on 
the Very Large Telescope at high spatial resolution (50-125 mas). With spatially re- 
solved 2 jjLTn spectra, we find that the H2 rotational temperature within the cometary 
knots is uniform. The rotational-vibrational temperature of the cometary knot (situ- 
ated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 
1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the 
central star) of the nebula (900 K), showing that the excitation temperature varies 
across the nebula. The obtained intensities are reasonably well fitted with 27 kms~^ 
C-type shock model. This ambient gas velocity is slightly higher than the observed 
[Hell] wind velocity of 13 kms~^. The gas excitation can also be reproduced with a 
PDR (photo dominant region) model, but this requires an order of magnitude higher 
UV radiation. Both models have limitations, highlighting the need for models that 
treats both hydrodynamical physics and the PDR. 

Key words: (ISM:) planetary nebulae: individual: NGC 7293 - ISM: jets and outflows 
- ISM: molecules - (stars:) circumstellar matter - ISM: clouds - Infrared: stars — 



1 INTRODUCTION 

In recent years it has become clear that knots of dense 
material are com mon in nebulae, in cluding Planetary 
Nebulae (PNe; e.g. lO'Dell et al.ll2002f) and star-forming 



* Based on observations with European Southern Observatory, 
Very Large Telescope with an instrument, SINFONI (the proposal 
numbers: 076.D-0807). 



regions (e.g. iMcCaughrean &: Mac Lowl I1997D . In the 
one of the best studied PN case, that of the Helix 
Nebula (NGC 7293) , knots in the inne r regions have 
a comet-like shape (jO'Dell &: HandronI il996) and are 
thus known as cometary knots. The Helix nebula is es- 
timate d to contain more th an 20,000 cometary-shaped 
knots (jMeixner et al.|[2005t ). The apparent commonal- 
ity of occurrence in PNe of these knots has lead to 
the assertion that all circumstellar nebulae are clumpy 
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in structure (e.s 



O'Dell et alJ 120021: ISoeck et al.ll2002t 



iMatsuura et al. l2005bD . As knots occupy the bright 



est parts of the planetary nebulae, understanding their 
physical nature is essential to understanding the domi- 
nant physics governing the nebula. 

The origin of these knots remains unknown. When 
the knots form is also disputed: they may have been 
present during the p receding AGB (As ymptotic Giant 
Branch) phase (e.g. iDvson et al.l Il989f ). or only have 
formed during the PN phase. The suggested formation 
mechanisms fall into two main scenarios. Firstly, photo- 
evaporation from condensations which pre -exist in the 
circum stellar envelope were considered. ISpeck et al.l 
(|2002f ) suggest that the UV radiation from the cen- 
tral star heats the surface of a condensation, emitting 
Ha and H2 on the facing side. Secondly, the interac- 
tion of a fast stellar wind with the slowly expanding 
AGB wind, during the early sta^e of PN phase, has 
been considered (e.g Vishniac 1994: 'Pittard et al.l l2005t 
iGarcfa-Segura et al.l r2006). Pittard ct al. (2005) found 



that cometary knots may be formed through instabili- 
ties where a supersonic wind impacts a subsonic wind. 

In addition to being seen as inhomogeneities within 
the ionized gas emissio n from nebulae, the kno ts also 
contain molecular gas (jSpeck et al. l l2002l l2003f). This 
yields a potential clue to their origin. iRedman et al.l 
(|2003|) argued that some molecular species formed in the 
AGB atmosphere can survive if molecules are shielded 
in clumps from UV radiations by extinction. But if 
the knots (re-)form during the ionized PN phase, the 
earlier AGB molecules were photo-dissociated, and the 
molecules reform within the knots, under conditions of 
relatively low density and higher UV intensity com- 
pared to the AGB wind. In these cond itions only sim ple 
molecules are expected to be present ([Woods! l2004l ). 

Using seeing- limited observations (1. 2 arcsec) of the 
H2 v= l-0 S(l) Hne in a single knot, iHuggins et al.l 



(|200 2l) found II2 gas in both the head and the tail of 
the knot, following the distribution of ionised gas traced 
by Ha and [Nil]. In contrast, CO thermal emission was 
found to emanate only from the tail. Higher angular 
resolu tion images obtained with HST (jMeixner et al.l 
|2005( ) resolved the globular H2 emission into crescent- 
shaped regions. With these previous studies, it is not 
clear whether the hydrogen molecules and ionised gas 
are co-located within the head or spatially separated. 

The excitation mechanisms for molecular hydro- 
gen in planetary nebula have long been contro- 
versial, with contention between fluorescent/thermal 
excitation in ph oton-dominant regions (PDR ) or 
shock e xcitation. Tielens fc HoU enbachl (|l985f ) and 
iBlack fc van Dishocck (1987) showed that vibrationally 
excited H2 (by FUV pumping) traces the surface 
of dense regions in the process of becoming photo- 
ionised. On the othe r hand, molecules can fo rm in 
a post-shock region (jNeufeld & Dalgarno 1989|) , and 
H2 can be excited by sho cks (jBeckwith ct ajj Il980l : 
iHollenbach fc McKed fl989[ ). The cooling region after 
the shock can be resolved if it is C-type (continuous) 
shock, but the region is too small to be resolved for J- 
type (jump) shocks. A comparison of models with spa- 



tially resolved spectra, which covers multiple line ratios, 
will help understanding the excitation mechanism of H2. 

The Helix Nebula (NGC 7293) is one of the nearest 
PNe, with a largest diameter of more than half a de- 
ree (iHora et al. 2006) and a parallax distance of 219 pc 
Harris et al.ll2007.) . Because of its distance, small-scale 
structures inside the nebula are well resolved and this 
nebula is used as a proxy to understand the structures 
found in PNe. 

We have observed a cometary knot in the He- 
lix Nebula, using the adaptive-optics-assisted, near- 
infrared integral field spectrometer on the Very Large 
Telescope. We have achieved 50-100 mas spatial res- 
olution. Our data provide the first spatially resolved 
spectra within a knot at 2 ^m and further this is the 
highest spatial resolution image-|-spectra of this PN at 
this wavelength. Our observations of H2 spectral line 
ratios allow us to derive the excitation temperatures 
within the knot, as well as differentiation between the 
possible excitation mechanisms of H2 . 



2 OBSERVATIONS AND ANALYSIS 

A knot in the Helix Nebula was observed by the Spec- 
trograph for INteg ral Field Observations in the Near In- 
frare d (SINFONI: [Eisenhauer et~al1l2003l: [Bonnet et al.l 
|2004 ) installed at the Cassegrain focus of the Very Large 
Telescope. The grating for the K-band was used. The 
wavelength band of the grating was mapped to the 2048 
pixel detector in the dispersion direction. We used two 
plate-scales: the spatial resolutions are 125x250 mas^ 
and 50x100 mas^. For all of the observations, a nearby 
star (RA=22:29:33.017 Dec=-20:48:12.68) was used as 
Adaptive Optics (AO) guide star. The spatial resolu- 
tion is usually determined by the pixel scale while us- 
ing the 125x250 mas^ plate-scale (i.e. the images are 
undersampled) , while for the 50x100 mas^ scale both 
the corrections of the AO system and the pixel scale 
are important. The spectral resolutions A/AA are 4490 
and 5090 for the 125x250 mas^ and 50x100 mas^ cam- 
eras, respectively. The observing log is summarized in 
Table [TJ The coordinates of our target knot Kl are 
22:29:33.414, -20:48:04.73 (J20 00) measured from the 
data us ed bvlO'Dell et all (l2004i) . This knot is noted as 
ID 6 in iMeaburn et al.l ([19981 ) and is about 2.5arcmin 
away from the central star (Fig.[T]). It is located at the 
inner rim of the ring-shaped nebula filled with cometary 
knots. 

The observing run was carried out on the first-half 
nights from 2"^ to 4'^ of November 2005. The sky was 
clear on Nov 3, with occasional thin cloud passing on 
November 2 and 4. The telluric standards were observed 
immediately after the target. The telluric standard stars 
were used for flux and spectral-response calibration, as 
well as for the measuren ients of the p oint spread func- 
tion. Model spectra from'P icklesI ([1998 1) are used as tem- 
plate spectra for the telluric standard stars. The wave- 
length resolution of the model spectra is 4000 only, but 
we assume that this difference in wavelength resolution 
is not critical for our analysis of the emission lines. 

The field of view was slightly jittered to minimize 
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Table 1. Observing log 



Target 


Pixel scale of the camera 
mas^ 


Pixel scale in reduced data 
mas^ 


Expt 


Telluric Standard 
Name Spec 


K-mag-f 


Kl 


125x250 


125x125 


300s X 10 


Hip 115329 


G2V 


7.552dz0.017 








600s X 10 










50x100 


50x50 


600s X 10 


Hip 023422 


G2V 


7.862±0.017 



t; Exp: Exposure time. 

t; K-mag: K'-band magnitudes from 2MASS 



the influence of bad pixels. The sky level was measured 
at +10 arcmin away in the declination direction, which 
is outside of the bright part of the nebula. Some residual 
of the sky level was left, and we use the integral field to 
remove the residual sky level. 

The ESO data reduction pipelin e for SINFONI 
on GASGANO (iModidiani et ai.ll200l was used. Dis- 
tortion correction and flat correction were adopted. 
Wavelength calibration used the Ne and Ar wave- 
length lamps, linearly interpolated throughout the en- 
tire wavelength coverage. After the data reduction with 
GASGANO, the final pixel was re-sampled to a 125 and 
50 mas scale for 125x250 mas^ and 50x100 mas^ plate 
scale images, respectively. 

The absolute flux calibration has some uncertainty, 
due to the possibility of occasional thin cirrus (up to 
0.2 mag) and the generic difflculty in the calibration of 
diffuse intensity using an AO-assisted point source. We 
reduced two 125 mas target spectra using two different 
exposure time independently, and the final spectra and 
images were obtained by average of these two data sets. 
The average of our measured flux over 2x2arcse c~^ is 



1.2 X 10 ergs cm sr 
the intensity of Kl at H2 



ISpeck et all (|2002f ) gives 
1-0 S(l) as - 5-10 X 



10~^ ergs""'^ cm~^ sr~^; the knots are not fully resolved 
(2 arcsec pixel" ^) and the flux of this knot is close to the 
detection limit 1 x 10~''^ergs~^ cm^'^ sr~^). A knot 
in the outer regio ns of the He l ix, w ith a similar flux 
level to our knot in Speck et "all (l200^ 's measurements, 
was observed bv lMeixner et al.1 (j2005D . They find a sur- 
face brightness of H2 of ~ 1 x 10""* ergs~^ cm~^ sr^^. 
We conservatively adopt a systematic calibration uncer- 
tainty of 50 percent. Further, the sky condition affects 
the correction of the atmospheric transmittance, espe- 
cially shortward of 2 /im and longward of 2.4 /im, due to 
the variation of terrestrial H2O. The relative intensities 
of lines below 2.0 fim and above 2.4 ^m have an error of 
~30%. We ignore the extinction at 2/im; as discussed 
later the extinction affects the absolute intensity less 
than 5% (^a=2. 128^111 = 0.05 mag). 



3 DESCRIPTION OF THE DATA 
3.1 Morphology 

Fig.[2]shows the image of the cometary knot Kl as seen 
in the three strongest H2 lines. The knot shows an clon- 



Figure 1. The location of th e cometary knot K l is plotted on 
the F658N ([Nn]+Ha) image dO'Dell et al."2004''). together with 
the AO guide star. The star symbol indicates the place of the 
central star. 



gated head with a narrower tail. The brightest emission 
is found in a crescent near the tip of the head. The 
crescent ends in two linear segments, indicated by the 
dotted lines in Fig. Hon the 2.12 /xm v=l-0 S(l) im- 
age. These segments are not co-aligned and deviate from 
the direction of the narrower tail. Overall, this gives the 
impression of a 'tadpole' shape, as opposed to the cylin- 
drical shape favoured by,0'Del l fc Handro n (1996), al- 
though in either case the overall shape is largely ax- 
isymmetric. Both the linear segments and the tail are 
brighter on the eastern side of the knot. 

The peak emission is located slightly behind the tip 
(Fig. 3]). There is a faint nebulosity around the bright 
head. This faint nebulosity appears not to be the wing 
of the point spread function (PSF). The faint halo ex- 
tends along the linear segments where the peak emission 
is much fainter. The higher resolution image (50 mas per 
pixel) in Fig. [3] also shows this faint nebulosity, suggest- 
ing that at least part of this faint nebulosity is real. 

The apparent diameter at the head is about 2.5 arc- 
sec including the faint rim, while the diameter decreases 
to about 2 arcsec along the tail. The transition from the 
linear segment to the narrower tail is visible as a kink, 
2.8 arcsec from the head on east side. It is less clearly 
visible on the west side. 
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Figure 2. Images of the cometary knot Kl, taken by 250 mas camera. The PSF at 2.12 /^m is found in (d). The direction of the central 
star is marked in (c) 



3.1.1 H2 emission from the surface of knots 

The deconvolved image (Fig. ^ shows a Hmb- 
brightened head of the knot. Within the head there may 
be an H2 empty region. FigiT] shows the radial cross- 
section of the head 1 arcsec inside from the brightest 
point. The intensity at the mid-point is about 65% of 
the peak in the raw image and about 40% in the decon- 
volved image. To demonstrate the shape, we modelled 
the radial profile in two dimensions, assuming a thin 
ring with a large hollow area inside (Fig[7| b and c). 
The diameter of the circular ring is 1.5 arcsec. The den- 
sity structure of the ring is assumed to be a Gaussian 
with a width of 0.075 arcsec. The radial cut in the de- 



convolved image is reasonably well reproduced by this 
model (Fig. [7] top). The H2 emitting region is probably 
a very thin surface of the knot, except at the tip. 



3.1.2 Comparisons with HST [Nil] and [OIII] images 

Fig. [5] shows the image of the knot Kl at 2.12 /im H2 
v=l-0 S(l), with a comparison of HST F658N (mainly 
[Nil] and some contribution of Ha) and F502N ([Olll]) 
images from O'Dell et al. (2004). The precise alignment 
between SINFONI and HST images is unknown, and the 
images were registered such that the tip of the knot is 
located at the same place at [Nil] and H2 v=l~-0 S(l). 
The [Olll] image shows the knot in absorption; the 



(a) L95 micron v= 1-0 S(3) 
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(b) 2. 1 2 micron v=l-0 S(l) 
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Figure 3. Image of the cometary knot Kl, taken with the 100 mas camera. The PSF at 2.12 fj.m is found in (c). 



tip is slightly off the peak fr om [Nil] image (F ip 
This has also been found by lO'Dell fc HennevI (I200C 
for other knots. 



A cross section of the H2 v=l-0 S(l) and [Nll] + 
Ha images shows that the decay of the intensities to- 
wards the tail is very fast (almost immediate) for the 
ionized lines, and slower for the H2 v=l-0 S(l). The 
main ionised region is a thin layer at the tip; the ma- 
jority of the material remains molecular or neutral. 



3.2 Spectra 

Fig. [TO] shows the spectra of the knot Kl at twelve 
regions indicated in Fig. [TTJ Within each region of 5 
by 5 pixels (except region (a) which has 5 by 4 pix- 
els), the SINFONI data have been averaged. At most 
twelve H2 lines (three of them with low signal-to-noise 
ratio) are detected by SINFONI. The intensities are 
summarised later in Table [H The strongest detected 
line is the 1.95 /xm H2 v=l-0 S(3), and the 2.12 fim 
H2 v=l-0 S(l) and 2.4 /im v=l-0 Q-branches are also 
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Figure 10. Spectra within tlie knot Kl at twelve 0.675x0.5 arcsec'^ areas defined in Fig. 1111 The y-axis shows the spectral intensity 
from —0.5 X 10** to 6 X 10** Jysr^^. The identifications of the H2 transitions are indicated in panel (a). 




Distance [arcsec] 

Figure 4. Annotated 2.12 image of the knot Kl (see text for 
details). 




n m, - 

1 2 J 4 5 

Figure 5. Deconvolved image at 2.12 ^m. A fain emitting region 
at the rim near the crescent is found, as well as narrowing of the 
tail. 
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Figure 12. The energy diagrams for area (a)— (1) defined in Fig. 1111 the column density for a given J upper level divided by the 
statistical weight as a function of upper state energy in cm~^. The square shows the measured energy levels. The thick lines show the 
LTE case for 1800 K and for column d ensities as given in the diagrams. For a comparison, we plot the 900 K LTE line, as derived for 
knots in the outer region of the nebula lICox et al llSii) in (a). The column density for the 900 K line is 8 X 10 cm . This temperature 
cannot fit the H2 line intensities measured in the knot Kl. 



strong. The H2 lines are unresolved at the resolution of 
R=5090. The 2.07 v=2-l S(3) and 2.15 Atm v=2-l 
S(2) 2.20 /im v=3-2 S(3) lines are only marginally de- 
tected. 

Several high transition lines within this wavelength 
range are not detected, such as the 2.00 /im v=2--l S(4), 
2.38 /im v=3-2 S(l). The Br7 line is not detected. 



3.2.1 Energy diagram 

Fig. [T2]shows the energy diagram of the H2 lines for sev- 
eral regions in the knot indicated in Fig. [TTJ Einstein co- 
efficients from I Turner et alj (|l977l ) are used. The slopes 
in the diagrams show that the line intensities follow a 
T = 1800 K Local Thermodynamic Equilibrium (LTE) 
distribution. At the bright area (a-c), the observed line 
intensities follow the 1800 K LTE line up to the upper 
energy of 8000 cm~^. The systematic uncertainty in the 
absolute H2 intensity (up to 50 per cent) affects the col- 
umn density, but not the excitation temperature which 
is determined by the slope. 



3.3 Line ratio map 



Fig. [13] shows the line ratio maps within the knot. 
Fig. [13] (a) maps the rotational temperature variations 
within the knot, represented by the 2.03 /im v=l-0 S(2) 
and 2.12 /im v=l-0 S(l) lines. This rotational temper- 
ature is uniform within the errors throughout the knot. 
This uniform rotational temperature is expected from 
the energy diagram (Fig. [12]), because all of the line 
intensities follow 1800 K LTE in any area. 



The vibrational temperature seems to vary within 
the knot but this is within the uncertainty. Fig. [13] (b) 
shows the line ratio of 2.24 /im H2 v=2-l S(l) and 
2.12 /im H2 v=l-0 S(l). This is the only combination to 
obtain the vibrational temperature within our observed 
data. The hue ratio tends to be highest at the tip (~ 
0.1) and to decrease towards the tail, although a careful 
treatment of v=l-0 S(0) line is needed. 
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Table 2. Intensities of H2 lines at area (a)— (c) and ratios with respect to the v=2— 1 S(l) line. 



Wav 


Transition 


Upper state 


energy 


Statistical weight 


Intensity /^xlO'^ 


Err§ 


Obs ratio^ 












(a) 


(b) 


(c) 




(a) 


(b) 


(c) 






K 






Wm 


^sr 


-1 










l.yoo 


1 n Qf'Q\t 

V — 1— u o^^oj ' 


ooDO 


Oolo 


00 


4.10 


2.35 


2.96 


U.lU 


220 


226 


223 


2.004 


v=2~l S(4) 


14764 


10262 


13 


<0.40§ 














2.034 


v=l-0 S(2) 


7584 


5271 


9 


0.67 


0.37 


0.49 


0.04 


36 


35 


36 


2.073 


v=2-l S(3) 


13890 


9654 


33 


0.08* 






0.09 


4 






2.128 


v=l-0 S(l) 


6956 


4834 


21 


1.86 


1.04 


1.32 


0.07 


100 


100 


100 


2.154 


v=2-l S(2) 


13150 


9139 


9 


0.07* 






0.08 


3 






2.201 


v=3-2 S(3) 


19086 


13265 


33 


0.02* 






0.03 








2.224 


v=l-0 S(0) 


6471 


4497 


5 


0.42 


0.24 


0.31 


0.06 


22 


23 


23 


2.248 


v=2-l S(l) 


12550 


8722 


21 


0.18 


0.09 


0.12 


0.04 


9 


9 


9 


2.386 


v=3-2 S(l) 


17818 


12384 


21 


<0.09§ 














2.407 


v=l-0 Q(l)t 


6149 


4273 


9 


1.84 


1.06 


1.30 


0.23 


99 


102 


97 


2.413 


v=l-0 Q(2)t 


6471 


4497 


5 


0.55 


0.29 


0.39 


0.12 


29 


28 


29 


2.424 


v=l-0 Q(3)t 


6956 


4834 


21 


1.85 


1.03 


1.27 


0.14 


99 


99 


96 


2.438 


v=l-0 Q(4)t 


7586 


5272 


9 


0.57 


0.31 


0.38 


0.18 


30 


30 


28 



^ Line ratio with respect to 2.128 /^m v=l— S(l) line 
^ 3(T of the noise level of the intensity in the error region indicated in Fig. 1111 
* These lines have a large (~30%) error in relative intensity calibration due to the influence of H2O in the terrestrial atmosphere. They 

are excluded from the comparison with models. 
* Marginal detections 
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Figure 13. The line ratio maps of the cometary knot Kl. The contour shows the surface brightness of 2.12 /xm H2 v=l-0 S(l) line at 
0.2 X 10-^, 0.5 X lO"'^, 1 X lO"'^, 2 x lO"'^ Wm-^sr"!. 



4 DISCUSSION 

4.1 Excitation diagrams and temperatures 

The H2 excitation diagrams are fitted by a single LTE 
temperatu r e of 1800 K. This is much higher than 
ICox et al.l (|l998t) who obtained 900 K by fitting pure- 
rotational li nes up to the S(7 ) transition. The region 
observed by ICox et ahl (jl99^ is located at the west- 
ern rim of the nebula, and is 5-6 arcmin away from the 



central star. The pixel size of their ISOCAM data was 
6 arcsec and emission from multiple knots contributed 
to each single pixel. In contrast, our target is a single, 
isolated knot located at the innermost region of the flock 
of knots (2.5 arcmin from the central star). 

lO'Dell et all (|2007[ ) also obtained an excitation 
temperature of 988 K from 5-15 /j , m spec tra in the outer 
part of the nebula. lO'Dell et all (|2007[ ) state that the 
distance of that slit positions is similar to that of 
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-4, I I I 

Figure 6. Cross cut of Fig. [2]along the radial axis, whose inten- 
sity is scaled to the maximum intensity of the knot. Approximate 
locations of cuts are indicated on the left panel. The data are 
smoothed with 250 mas bin along long axis of the knot. The in- 
tensity peaks on east side (right side in this figure) of the knot 
are marked with stars in each panel. The bottom panel shows 
the cross cut of the telluric standard (bold line), airy disk for 8.2 
meter circular aperture (thin solid line), and pixel size (125 and 
250 mas; thin dash lines). The spatial resolution is determined 
by the pixel scale (sampling rate). The star marks traces that the 
head is widen until at distance {d)=2.0 arcsec and then narrowing 
follows d >2.5. Tail of the PSF is less than 10% of the peak at 
0.5 arcsec. 



I Cox et all (|l998t) . Although they do not provide de- 
tailed information on the slit positions and observing 
mode for their two Spitzer spectra, the Spitzer archive 
can provide this, suggesting their two slit positions are 
5.6 and 4.2arcmin from the central star, and the slit 
size is 3.6x57 arcsec^. Further three spectra at different 
slit positions are available in arc hive, but th ese do not 
fit the description of the data in lO'Dell erai.i (j2007f ). 

The difference in measured temperatures indicates 
that the excitation temperature of H2 is not uni- 
form within the Helix nebula: the H2 molecules reach 
higher temperature within the inner region. From about 
2.5arcmin to 5arcmin away from the central star, the 
excitation temperature decreases from 1800 K to 900- 
1000 K. This provides evidence for temperature varia- 
tions within the nebula. 
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Figure 7. (a) Cross cut of the observed intensity and at 2.12 /im 
and PSF deconvolved profile along the radial direction in the 
head. The intensities are normalised at the peak. The location is 
1 arcsec away from the brightest tip. Model profile in panel (a) 
is to demonstrate hollow shape in projected image if the knot is 
seen from the central star (panel b) with input parameters of the 
shape (panel c) 



iHora et all (|2006D show that the [3.6] -[4.5] vs 
[4.5] — [8.0] colour varies within the nebula. All three 
bands are dominated by H2 lines. Knots located at the 
inner rim of the main nebula are brighter in the 4.5 /xm 
band. The strongest expected line in this band is the 
4.69 /im H2 0-0 S(9), whose upper energy is the high- 
est among the dominant H2 lines w ithin the IRAC fil- 
ters. The colour-colour diagram of IHora et all (gooi) 
also suggests a higher excitation temperature in the in- 
ner region, and global variation in the excitation tem- 
perature in the entire nebula. Analysis of spectra with 
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Figure 8. Comparison of tlie SINFONI 2.12 fim image and HST 
images. The SINFONI and HST images were aligned on the tip 
of the knot. Log intensity scale is used for all of the images. The 
size of the image is 22.3 X S.Sarcsec'^ 



several slit positions within a nebula shows that the H2 
temperature is not un iform within a PN (|Hora et al.l 
I1999I : [Pavis et al.ll2003D . Excitation temperature varia- 
tions appear normal within a single PN. 

Furthermore, H2 excitation temperatures in PNe 
have been found as high as 2000 K (Hora et 
al. 19 99; Davis et al. 2003) and lower than 
lOOOK dCox et aL'^jgi: ' Bernard-Salas fc TielensI [20051 : 
iMatsuura fc Ziilstra 2005^ Variations in excitation 
temperature both between PNe and within a single PN 
appear to be normal 



4.2 Column density 

The H2 line intensities are fitted using column densi- 
ties in the range 1 x 10^^-8 x lO^^ cm"^. This range 
is smaller than the ICox et al] (|l998l ) measurements of 
3 X 10^* cm^^. Our target is an isolated knot, while the 
low spatial resolution of the ISOCAM data almost in- 
cluded multiple knots within the field of view, h ence the 
H2 intensity is much higher (jSpeck et al.ll200^ . 

The column density of H2 can be converted to a hy- 
drogen mass of 2 X 10~®Mq. Here we assume the knot 
has a column density of 4 x 10^^ cm~^ (i.e. to emit aver- 
age intensity over 2x2 arcsec^ at the excitation temper- 
ature of 1800 K), the distance is 219 pc, the dimension 
of the knot is 2x2 arcsec^ and the density is uniform 
within this area. The estimated hyd rogen mass is a fac- 
tor of 1000 less than the estimate of lO'Dell fc HandronI 
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Figure 9. Cross cut of Fig. |8]along the long axis. 



(199(f), who used the extinction of the [Oil] line as a 
mass tracer. Our infrared H2 lines trace only highly 
excited H2, as expected from the upper state energy 
(>60 00K). Co lder H2 gas is inefficient at emitting these 
lines (|Cox et a l. 1998). 

In order to measure the mass of cold H2 gas more di- 
rectly, we need observations at UV wavelengths, where 
H2 lines could be found in absorption. However, this 
approach may be compromised by the contribution to 
H 2 absorption by the ISM Adopting the hydrogen mass 
of lO'Dell fc HandronI (|l996l ) with a comparison of our 
estimated H2 mass indicates that H2 is heated only at 
the surface, but is in fact present throughout the knot. 
This favours a scenario where the detected H2 was al- 
ready present within the knot, i.e., it is not necessary 
to assume that it formed from chemical reactions from 
atomic/neutral hydrogen within this surface. However, 
this argument is based on an expected correlation be- 
tween dust extinction and H2. A direct detection is re- 
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Figure 11. The twelve regions shown by boxes on tlic 2.12 fim 
image. The spectra are averaged within that area, which corre- 
sponds to the diagram in Fig. 1121 The right top box is the area 
used to estimate the error of the intensities. 



quired to confirm the presence of a large, cold H2 reser- 



4.3 Excitation mechanisms of molecular hydrogen 

The cometary knot Kl emits a 1-0 S(l) line intensity 
of 2 X 10~^ Wm~^ sterad"^. The adopted uncertainty 
is a factor of two. The 2-1/1-0 S(l) line ratio is ~ 0.1 ± 
0.02 at the tip, using the uncertainty on the response 

cor rection. 

iBurton et al.l ()1992[ ) calculated the line intensities 
of 1-0 S(l) and 2-1 S(l) hues for J-type shocks, C- 
type shocks and in PDR. We found that within their 
model the measured line intensities and ratios can be 
fit equally well under several model conditions. 

• A photo-dissociation region (PDR) heated by UV 
radiation where the density is n = 10^ cm""^ and the 
UV strength is Go = 1.2 x 10^ (within the range of 
Go =1-1.2x10^) 

• C-type shocks with upstream density uq = 
Iff' cm^"^ and shock velocity Vg ~ 27kms~' (within the 
range of 27-28 km s~^). 

• J-type shocks no = 10® cm^'^ and shock velocity 
Vs - 9kms-H9-10kms"i). 

Although there is some uncertainty in absolute in- 
tensity, the line ratio is the strongest constraint on the 
abo ve conditions. 

iMeaburn et al.l ()1998f ) estimated the molecular hy- 
drogen density of another knot to be 8.9 x 10^ cm^'^ 
using the extinction in the [Olll] 5007 A line. This knot 
appears brighter than Kl, and larger in radius when 
observed in the [Nil] line. We follow their method to 
estimate the density of the knot Kl. The absorption at 
[Olll] is 0.73 with respect to the continuum fSect l3.lT^ . 
This corresponds to an extinction coefficient c = 0.15, 
E{B — V) =0.14 mag, and an H2 column density of 
4. Ox 10^° cm^^. If the diameter of the knot is 1.5 arcsec, 
an H2 number density of 8 x 10^ cm~^ is obtained for a 



distance of 219 pc. The derived H2 density is consisten t 
for PDR model of H2 excitation in lBurton et al.l (|l992t ). 

The FUV (6-13.6 eV) flux at the knot Kl is about 
Go = 8, wher e Go is the F UV radiation measured 
in units of the iHabingj (jl968D flux. This is based on 
ISu et al.l ()2007D 's estimate assuming a luminosity of the 
central star of 76 Lq ai id a distance of 219 pc. The re- 
quired UV radiation in iBurton et a l.' (1992) 's model is 
Go > 1 X 10^, which is an order of magnitude higher 
than the estimated UV radiation field stren gth at the 
knot Kl. We used the UCL-PDR model jBell et al.l 
I2OO5) to calculate the conditions independently. At an 
UV radiation fiel d of Go = 8, we obtain a temperature 
below 100 K. RoUig et al.l (|2007D compare benchmark 
calculations of indepen dent PDR codes, including t he 
UCL-PDR model and iTielens & HoUenbach (l^'s 
model that was incorporated in Burton ct al. (1992)'s 
H2 model. They find consistent gas temperatures among 
PDR models in their benchmark calculations. These re- 
sults suggest that the PDR models can reproduce the 
measured LTE gas temperature (1800 K) of the H2 lines 
in the knot Kl only, if the UV radiation stre ngth is in- 
crease d. A similar conc l usion is derived by ICox et ah! 
(|T99l and 10 'Dell et"all (pOOl . The value for Go is de- 
fined as the flux between 6-13.6 A a narrow-band UV 
radiation field, adopted for interstellar radiation field. 
If we consider the continuous stellar spectrum of the 
central star, the U V flux increases by a factor of 10 
(jO'Dell et al.ll2007D . However, this remains a factor of 
250 short of the required flux. 

The speed of ambient gas for C-type shocks is 
Vg ~ 27kms~' to fit the II2 lin e intensities and lin e 
ratios, from theoretical work by IBurton et al.l (|1992f ). 
The required velocity can vary by ~ ilOkms"-'^ de- 
pending on the assumed magnetic field strength and 
the iron fraction. The [Hell] and [Olll] lines show 
an expa nsion velocity of 13kms~' near the central 
star (Me aburn et al.ir2005t ). Within an expanding neb- 
ula, hydrodynamic eff ects will c ause significant velocity 
gradient s (iSchocnberner et al.l [2005). as are observed 
in PNe (|Gesicki et alT 20031 ). The overpressure of the 
ionised region dominates, and (if present) the pressure 
from the inner, hot bubble might be added. The same 
processes wi ll occur for e ach knot at the ionized, fac- 
ing edge. iMeaburn et al.l (| 19981 12005[ ) postulate turbu- 
lent velocities ^ lOkms"' (i.e. larger than the sound 
speed). A C-type shock velocity could be associated 
with such motion. The gas density of PNe is typically 
uq ~ 10'*cm~^. This is consistent with the required 
upstream density. C-type shock excitation of H2 is pos- 
sible. Stronger observational constraints on the velocity 
structure of the knot would be helpful. 

The presence of magnetic fields has been reported 
in several you ng or pre-planetar y nebulae using sub- mm 
polarization (jSabin et al.ll2007l ). They find a magnetic 
field of ~lmG at 5 x 10^® cm from the central star in 
other PNe. The Helix knots are located 10 times further 
from the star than their measured location. A dipole 
field would decay as r~^: this would leave a negligible 
field in the Helix. A solar-type field decays as r~^, and 
a frozen toroidal field (as favoured bv lSabin et al.|[2007l ) 
may decay slower with radius. The formation of a knot 
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may strengthen after its embedded magnetic field. The 
field required for the C-type shock is plausible within 
the knot, compared to the stronger fields detected in 
more compact shells. The inter-knot medium is likely 
to show a weaker field. 

An upstream (wind) density of 10® cm~^ is re- 
quired for the J-type shock model. Although the den- 
sity with in knots themselves are recorded at lO^cm"'^ 
(jO'Dell fc Handronl99& : ,Meaburn et al.lll998l ). it is un- 
likely that the upstream region has such a high density. 



5 CONCLUSIONS 

We have investigated the detailed structure of a single 
knot close to the inner edge of the main ring of the He- 
lix nebula. We find that the rotational- vibrational tem- 
perature of H2 is as high as 1800 K for this innermost 
cometary knot. The rotational temperature is uniform 
within the knot, and the vibrational temperature ap- 
pears to follow the same distribution except for a possi- 
ble decrease towards the tail. The derived temperature 
is much higher than previously measured (900 K) for 
knots in the outer region. The excitation temperature 
changes with radial distance from the central star. The 
studied knot has a wide head, with the H2 distributed 
in a crescent. 

We examine the possible molecular hydrogen ex- 
citation mechanisms. Based on the line intensities and 
ratios, C-type shocks can provide a fit using plausible lo- 
cal conditions, althought the required velocity is slightly 
higher than the observed one. This excitation requires 
the presence of an as-yet undetected magnetic field. A 
J-type shock model does not fit. We also tested the pos- 
sibility that H2 lines are excited in PDR, but PDR mod- 
els have difficulties to reproduce the very high temper- 
ature, unless the UV field is increased by two orders of 
magnitude over the stellar radiation field. A further test 
of the models will require observations of the velocity 
structure of the head and tail. We also lack a synthe- 
sis model combining collisional and UV excitation, with 
hydrodynamical interactions. 
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